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. Abstract. The elliptic flow parameter for Kg and A + A has been measured at 

\ mid-rapidity in Au + Au collisions at y/s NN = 200 GeV by the STAR collaboration. 

The v 2 values for both K s and A + A saturate at moderate pr, deviating from the 
hydrodynamic behavior observed in the lower p? region. The saturated v 2 values 
and the px scales where the deviation begins are particle dependent. The particle- 
qq . type dependence of v 2 shows features expected from the hadronization of a partonic 

ellipsoid by coalescence of co-moving quarks. These results will be discussed in relation 
■ to the nuclear modification factor (Rcp) which has also been measured for K s and 

q ! A + A by the STAR collaboration. 
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1. Introduction 

The elliptic component of the event-wise azimuthal anisotropy of particle production 
(i.e. elliptic flow or v 2 ) is thought to probe the early stages of relativistic heavy ion 
collisions j^j- Measurements of v 2 for identified particles 0111 and charged hadrons EJ 
at the Relativistic Heavy Ion Collider (RHIC) indicate a conversion of spatial anisotropy 
to momentum anisotropy near the hydrodynamical limit IB] • For px greater than 
2 GeV/c, however, the charged hadron v 2 deviates from hydrodynamical calculations 
and saturates at a large value approximately independent of pr up to 6 GeV/c 6\. The 
measurements of v 2 for Kg and A + A also indicate a saturation at moderate px- Models 
using large parton energy loss [TDJ and transport opacity jTT] have been discussed in 
relation to the saturation and centrality dependence of charged hadron t> 2 at large px- A 
saturated v 2 could also arise if the extent of the pr region where soft processes dominate 
the spectrum is particle dependent [T2]. In this case the hyperon v 2 may continue to 
rise at intermediate px (perhaps following hydrodynamic model calculations) while v 2 
of the less massive meson decreases. It has also been suggested that if a partonic state 
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exists prior to hadronization, the process of particle formation at moderately high px, 
by string fragmentation, parton fragmentation ^3] or quark coalescence [THITo1ITHl ll7j. 
may lead to a dependence of v 2 and Raa on particle type. As such, it's possible that 
these measurements will provide information on the existence and nature of an early 
partonic state. 

In this Letter, we report the measurement of v 2 at mid-rapidity, \y\ < 1.0, for K s 
and lambda (A) + antilambda (A) with 0.4 < p T < 6.0 GeV/c and 0.6 < p T < 6.0 GeV/c 
respectively in Au+Au collisions at y/s NN = 200 GeV. 

2. Analysis and Results 

This analysis uses 1.6 x 10 6 minimum-bias trigger events and 1.5 x 10 6 central trigger 
events detected in the STAR detector system [TH]. The particles, K s , A and A, were 
identified from the charged daughter tracks produced in the decays K s — > n + + tt~, 
A — > p + 7r~ and A — > p + 7r + . A detailed description of the analysis, such as track 
finding, decay vertex topology cuts, and the estimation of detection efficiency, is given 
in Refs. jHHJ. 

We use the yield as a function of {4>%j — ^f) to calculate v 2 = (cos[2(0jj — ^f)}), 
where 0y is the azimuthal emission angle of particle i in event j and tyf is the reaction 
plane angle for event j, where, to remove autocorrelations, the decay daughter tracks 
associated with particle i are excluded from its calculation. Within statistical errors, 
the A v 2 is the same as A v 2 , so they are summed together. 



Table 1. The systematic errors for minimum bias vi from non-flow effects (n-f) and 
background contamination (bkg). The values represent the absolute errors. The px 
resolution, 5px/px is also listed. 







A + A 


p T (GeV/c) 


1.0 


2.5 


4.0 


1.0 


2.5 


4.0 


v 2 (bkg) 


+0.000 -0.001 


+0.001 -0.007 


+0.003 -0.018 


+0.001 -0.007 


+0.005 -0.001 


+0.005 -0.001 


v 2 (n-f) 


+0.00 -0.01 


+0.00 -0.04 


+0.00 -0.03 


+0.00 -0.01 


+0.00 -0.04 


+0.00 -0.04 


5p T /pT 


0.016 


0.027 


0.037 


0.016 


0.027 


0.037 



Possible sources of systematic error in the calculation of v 2 are correlations unrelated 
to the reaction plane (non-flow effects), uncertainties in the extraction of yields from 
the invariant mass distributions, the particle momentum resolution {Spt/pt), and biases 
introduced by the cuts used in the analysis. Tabled lists the dominant systematic errors 
for three transverse momenta. The non-flow systematic error is dominant. The non-flow 
effects for charged particle v 2 are discussed in Refs. jHl Ej but, the particle dependence 
of these effects has not been measured. We assume a similar magnitude of non-flow 
contribution to A + A and K s v 2 - A 4-particle cumulant analysis of A + A and K s v 2 
will be less sensitive to non-flow effects but, to be conclusive, will require a larger data 
sample than is currently available. 
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STAR Preliminary (Au+Au;>|5^= 200 GeV; Rapidity |y| < 1.0) 
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Figure 1. The elliptic flow of Kg and A + A as a function of p T for 30-70% (a), 
5-30% (b), 0-5% (c) and 0-80% (d) of the collision cross section. The error bars shown 
are statistical errors only. 

Fig. Q shows v 2 of Kg and A + A as a function of pt for the centrality intervals, 
30-70% (a), 5-30% (b), 0-5% (c), and 0-80% (d) of the geometrical cross section.The p T 
dependence of the v% for all the centrality bins has a similar trend. There is a saturation 
and particle dependence at moderate px for each of the centrality intervals and, as such, 
the saturation for the minimum-bias (0-80%) cannot be due to the superposition 
of drastically different pt dependencies in various centrality bins. This measurement 
establishes the saturation and particle dependence of v 2 at the moderate to high p T 
region. 



STAR Preliminary (Au+Aujji - = 200 GeV; |y|<1.0) 
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Figure 2. (Top panel) The elliptic flow of Kg and A + A as a function of mj — m 
for 30-70% (a), 5-30% (b) and 0-5% (c) of the collision cross section. The error 
bars shown are statistical errors only, (bottom panel) The same scaled by the initial 
eccentricity. 

Fig. Eland Fig. 01 show the same data as Fig. ^ but plotted versus rriT — mo which is 
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STAR Preliminary (Au+Au;^s~~ = 200 GeV; Mid-rapidity) 




Figure 3. The elliptic flow of K s and A + A as a function of px (top panel) and 
rriT — mo (bottom panel) for 0-80% of the collision cross section. The error bars shown 
are statistical errors only. 



approximately the kinetic energy of the hadron. In the low momentum region the K s and 
A+A t>2 appears to fall on a single straight line. The hydrodynamic calculations in Fig. El 
seem to capture this trend. In the bottom panels of Fig. El we've scaled the v 2 values by 
the eccentricity of the overlap region for the various centralities. Hydrodynamic models 
predict that v 2 should scale with this initial eccentricity. 

In Fig. |U (top) we show v 2 (j>t) for Kg, A + A, and charged hadrons [20] along with 
hydrodynamic model calculations of v 2 for identified particles [7|. Below pt ~ 1.2 GeV/c 
v 2 is consistent with the calculations and in agreement with the previous results for Kg 
and A + A at y/s NN = 130 GeV [4 . Contrary, however, to hydrodynamical calculations, 
where at a given p? heavier particles will have smaller v 2 values, the measured v 2 of 
the heavier hyperon saturates at a value significantly larger than the v 2 of the lighter 
K s meson. The px scale where the measured v 2 deviates from the hydrodynamical 
prediction is particle dependent with the hyperon v 2 following the prediction up to 
Pt ~ 2.0 GeV/c while the K s v 2 deviates much sooner. 

The ratio (Rcp) of the yields in central and peripheral collisions scaled by the 
number of binary nucleon-nucleon collisions (Nu nary ), may also be sensitive to the effects 
of energy loss and hadronization via parton coalescence ^21 E] • I n Fig. 0] (bottom) we 
show Rcp for K s , A+A and charged hadrons using the centrality intervals 0-5% (central) 
and 60-80% (peripheral) [2UI221- The charged hadron spectrum at pt > 2 GeV/c for 
the 60-80% centrality bin approximately follows binary collision scaling without medium 
modification |22j- As such, when this bin is used, Rcp approximates Raa- The bands 
in Fig. |U represent the expected values of Rcp for N binary and N part scaling including 
systematic variations within the calculation |22j . 

For p T < 5 GeV/c, the K s and A + A yields are suppressed (relative to N binary 
scaling) by different magnitudes. In addition, the pt scale associated with the onset 
suppression has a dependence on particle-type that is similar to the dependence in t> 2 
for the onset of saturation. At pt ~ 5.0 GeV/c, Rcp values for K s and A + A are both 
approaching the value of the charged hadron Rcp- 
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Figure 4. Minimum-bias t>2 (top panel) and Rcp (bottom panel) for Kg and A + A. 
The vi errors are statistical only. The widths of the gray bands represent uncertainties 
in the model calculation of Nbi nary and N part . Charged hadron V2 and Rcp are also 
shown. 

3. Discussion 

Although Rcp depends only on the yield in the central and peripheral bins, and the 
v 2 in Fig. |3] is from a minimum-bias centrality interval, the two parameters may be 
intimately related. The differential elliptic flow (^(pr)) measures the ratio of the pt 
spectrum of particles emitted in the direction of the reaction plane (in-plane) to that 
of particles emitted perpendicular the reaction plane (out-of-plane). In hydrodynamic 
models, we expect the pressure gradient to be larger in the in-plane direction than 
the out-of-plane direction. In this case, v-i will be the ratio of a pt spectrum with 
more hydrodynamic flow to one with less flow. As such, taking the ratio (Rcp) of a 
Pt spectrum that exhibits large flow (central) and one exhibiting less flow (peripheral) 
should lead to a very similar px and particle-type dependence. A detailed study of v-i 
and Rcp for identified particles should reveal the extent (in p T and centrality) to which 
hydrodynamical models are valid. 

A surface emission scenario, where partons traversing a dense medium experience 
large energy losses, has been discussed in relation to the large, pr independent t>2 
measured for charged hadrons [XOj - This mechanism should also lead to a suppression of 
particle production in central Au+Au collisions. This scenario, however, is inconsistent 
with STAR measurements of t>2 and Rcp for Kg and A + A. The smaller suppression 
manifested in the A + A Rcp contradicts the larger azimuthal anisotropies manifested 
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in A + A V2 values. In addition, calculations based on a surface emission model ^U] 
cannot produce v 2 values as large as those measured. 

The absence of a net suppression of A + A for p T from 1.8-3.5 GeV/c in central 
Au+Au collisions could also indicate the presence of dynamics beyond the framework 
of parton energy loss followed by fragmentation. The stronger dependence on centrality 
(and thus parton density) for baryon production indicated by the larger Rqp would 
naturally be expected from multi-parton mechanisms such as gluon junctions quark 
coalescence jTU, or recombination [T^]. Within the framework of these models, the 
measured v 2 and Rep features may reflect the anisotropy and hadronization properties 
of the bulk quark matter. Fig. shows v 2 of Kg and A + A as a function of pr where 
the V2 and px values have been scaled by the number of constituent quarks (n). Above 
p T /n ~ 0.8 GeV/c, the U2/11 vs pr/ n is the same, within errors, for both species. In a 
scenario where hadrons at intermediate pp (~ 1 — 5 GeV/c) are predominantly formed 
from bulk partonic matter by quark coalescence, e.g. Ref. |14j . t^/n should reveal the t>2 
developed by partons prior to the hadronic phase. The verification of this scenario, to 
the exclusion of other possible explanations, would be strong evidence for the formation 
of a quark-gluon plasma at RHIC. 



STAR Preliminary (Au+Au; \|s^ = 200 GeV; |y|<1 .0) 
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Figure 5. The v-i parameter for K% and A + A scaled by the number of constituent 
quarks (n) and plotted verses pxjn. 

The scenario discussed in Ref. fT2j . where soft processes dominate the hyperon 
t> 2 up to a higher p T than mesons, could also lead to a particle dependence for R C p 
that is qualitatively consistent with these measurements. Quantitative calculations of 
V2 and predictions for Rcp from this scenario, however, are still needed. Up-coming 
measurements of Rcp for identified particles in d + Au collisions at RHIC will also 
make it possible to study the effect of initial state interactions on Rcp- 
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4. Summary 

We have reported the measurement of i>2 for Pt up to ~ 6.0 GeV/c for ifij and A + A 
from Au + Au collisions at y/s NN = 200 GeV. For p? < 1.2 GeV/c, hydrodynamic model 
calculations agree well with the pt and mass dependence of the measured t>2. At this 
low momentum region K s and A + A v 2 lie on a single straight line when plotted versus 
— m . In the moderate pt region, however, the particle type and pt dependence 
of V2 suggests hydrodynamics no longer describes the collision dynamics. The value of 
V2 for Kg saturates earlier and at a lower value than the A + A V2- Measurements of 
Rcp show that the suppression of particle production in central collisions depends on 
particle-type in a similar way. The measurement of the particle-type and pt dependence 
of t> 2 and Rcp at moderate pt may provide a unique means to establish the existence 
(and study the properties) of a quark-gluon plasma that may be formed in collisions at 
RHIC. 

Acknowledgments: This work has been partially supported by NSF grant PHY- 
03-11859. 



5. References 



H. Caines (for the STAR Collaboration), these proceedings. 

J.-Y. Ollitrault, Phys. Rev. D 46, 229(1992); H. Sorge, Phys. Rev. Lett. 82, 2048(1999). 

C. Adler et al. (STAR Collaboration), Phys. Rev. Lett. 87, 182301 (2001); S. Esumi (for the 

PHENIX Collaboration), e-print jnucT-ex/0210012| 
C. Adler et al. (STAR Collaboration), Phys. Rev. Lett. 89, 132301 (2002); P. Sorensen (for the 

STAR Collaboration), J. Phys. G 28, 2089 (2002). 
C. Adler et al. (STAR Collaboration), Phys. Rev. C 66, 034904 (2002). 

C. Adler et al. (STAR Collaboration), Phys. Rev. Lett. 90, 032301 (2003). 

P. Huovinen, P.F. Kolb, U. Heinz, P.V. Ruuskanen, and S.A. Voloshin, Phys. Lett. B 503, 58(2001). 

D. Teaney, J. Lauret, and E.V. Shuryak, Phys. Rev. Lett. 86, 4783(2001). 

E. Wang and X. N. Wang, Phys. Rev. Lett. 89, 162301 (2002); 
E. V. Shuryak, Phys. Rev. C 66, 027902 (2002). 

D. Molnar and M. Gyulassy, Nucl. Phys. A 697, 495 (2002) [Erratum-ibid. A 703, 893 (2002)] 
M. Gyulassy, I. Vitev, X. N. Wang and P. Huovinen, Phys. Lett. B 526, 301 (2002). 
ZW. Lin and CM. Ko, Phys. Rev. C 65, 034904(2002). 
D. Molnar and S. A. Voloshin, e-print|nucUh/03020lJ 



Z. W. Lin and C. M. Ko, Phys. Rev. Lett. 89, 202302 (2002) 

R. J. Fries, B. Muller, C. Nonaka and S. A. Bass, e-print|n ucl-th/0301087 

V. Greco, C. M. Ko and P. Levai, e-print: nucl-th/0301093 



C. Adler et al. (STAR Collaboration), Nucl. Instr. Meth. A 499, 624 (2003). 
C. Adler et al. (STAR Collaboration), Phys. Rev. Lett. 89, 092301 (2002); K. Safarik, I. Kraus, 
J. Newby and P. Sorensen, AIP Conf. Proc. 631, 377 (2003). 



K. Filimonov (for the STAR Collaboration), e-print nucl-ex/0210027 
H. Long (for the STAR Collaboration), these proceedings. 
J. L. Klay (for the STAR Collaboration), e-print nucl-ex/0210026 
S.E. Vance and M. Gyulassy, Phys. Rev. Lett. 83, 1735 (1999). 



